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We compute the /iQs and pion form factors using partially twisted boundary conditions. The 
twists are chosen so that the K(^t, form factors are calculated directly at zero momentum transfer 
[q^ = 0), removing the need for a interpolation, while the pion form factor is determined at 
values of close to q^ = 0. The simulations are performed on an ensemble of the RBC/UKQCD 
collaboration's gauge configurations with Domain Wall Fermions and the Iwaski gauge action 
with an inverse lattice spacing of 1.73(3)GeV. Simulating at a single pion mass of 330 Me V, we 
find the pion charge radius to be (r )330MeV = 

0.354(3 l)fm2 which, using NLO SU(2) chiral 
perturbation theory, translates to a value of (r^) = 0.418(31)fm^ for a physical pion. For the 
value of the K(t, form factor, f^^{q^), determined directly at q'^ = 0, we find a value of f^j^ifi) — 
0.9742(41) at this particular quark mass, which agrees well with our earlier result (0.9774(35)) 
obtained using the standard, indirect method. 
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1. Introduction 
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IV^^I =0.97408(26) 

f^^(O) =0,9644(47) 
yf^=1. 198(10) 
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Over the last two years as part of our Domain Wall Fermion (DWF) physics programme we 
have been looking at the K —>■ n£V( (K^^) form factor at zero momentum transfer. Since the experi- 
mental rate for ^^^3 decays is proportional to |Kip|/^^(0)p, a lattice calculation of the form factor, 
fini^^) ~ 0' provides an excellent avenue for the determination of the Cabibbo-Kobayashi- 
Maskawa (CKM) [1] quark mixing matrix element, \Vus\- 

The uncertainty in the unitarity relation of the CKM matrix \Vud^ + \Vus^ = 1 (we ignore 
\'Vub\ since this is very small), is dominated by the precision of \Vus\- In Fig. |l| we show the latest 
determinations of \Vud\ [2] and \Vus \ [3]. For comparison, we also show the unitarity relation. Since 
it is important to establish unitarity with the best precision possible, it is essential that we decrease 
the error in \Vus\- 

The value of f^j^{0) used in determining \Vus\ 
in figure [I] was determined using standard meth- 
ods [4, 5] involving periodic boundary conditions 
in the recent paper [3]. There, the Kgi, form fac- 
tor is calculated at q^^^ = {rriK — rrijiY and sev- 
eral negative values of for a variety of quark 
masses. This allows for an interpolation of the 
results to q^ = 0. The form factor is then chi- 
rally extrapolated to the physical pion and kaon 
masses. The final result for /^^^(O) quoted is then 
[3] /^^(O) = 0.9644(33) (34) (14) where the first 
error is statistical, and the second and third are es- 
timates of the systematic errors due to the choice 
of parametrisation for the interpolation and lattice 
artefacts, respectively. This gives us a value of 
|K.v|= 0.2249(14). 

More recently, we have developed a method that uses partially twisted boundary conditions to 
calculate the K^t, form factor directly at 17^ = [6], thereby removing the systematic error due to 
the choice of parametrisation for the interpolation in q^. We have also used partially twisted bc's 
to calculate the pion form factor at values of q^ below the minimum value obtainable with periodic 
bc's. In contrast to recent studies this allows for a direct evaluation of the charge radius of the pion. 
The method was developed and tested in [6] and now applied in a simulation with parameters much 
closer to the physical point. 

In this paper we discuss our findings for the pion form factor from [7] and our progress in im- 
proving the precision of our result for f^„{0) from [3] using partially twisted boundary conditions. 
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Figure 1: Bands showing the current limits on 

\V„d\ [2], and \Vus\ [3]. 



2. Simulation Parameters 



The computations are performed using an ensemble with light quark mass am„ = am^ = 0.005 
and strange quark mass arris = 0.04 from a set of A// = 2 + 1 flavour DWF configurations with 
(L/a)3 xT/axLs = 24^ x 64 x 16 which were jointly generated by the UKQCD/RBC coUabora- 
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tions [8] using the QCDOC computer. The gauge configurations were generated with the Iwasaki 
gauge action with an inverse lattice spacing of a"' = 1.729(28)GeV. The resulting pion and kaon 
masses are m„ k, 330MeV and mx ~ 575MeV, respectively. 

In tliis work we use single time-slice stochastic sources [10], for which the elements of the 
source are randomly drawn from a distribution & = Z(2) ® Z(2) which contains random Z(2) 
numbers in both its real and imaginary parts. With sources of this form we find that the computa- 
tional cost of calculating quark propagators is reduced by a factor of 12. For more details on the 
simulations, see [7]. 

3. The Form Factors 

Here we briefly outline the main features of our method and we refer the reader to our earlier 
papers for more details [3,6,7]. 

The matrix element of the vector current between initial and final state pseudoscalar mesons 
Pi and Pf, is in general decomposed into two invariant form factors: 

(/V(/'/)I^Ml^K/'/))=//tp,(^')(A+/'/)M+//^^p/^')(A--/'/)M> (3-1) 

where = —Q^ = [jpi — pff-. For K ^ n, = sy^u, Pi = K and Pf = K. For tt ^ tt, = 
^uy^u — ^dy^d, Pi = Pf = n and from vector current conservation, fj^iq^) = 0. The form factors 
fp.p^{q^) and fp.p^{q^) contain the non-perturbative QCD effects and hence are ideally suited for a 
determination in lattice QCD. 

In a finite volume with spatial extent L and periodic boundary conditions for the quark fields, 
momenta are discretised in units of 2n/L. As a result, the minimum non-zero value of for the 
pion form factor in our simulation is k. —0.15 GeV^, while for the K ^ n form factor 

q^ = {EKm-En{Pf)? - {Pi-Pff ■ (3.2) 

For Pi = and 2n/L with pf = 0, we have q^ ?a 0.06 GeV^ and —0.05 GeV^, respectively, pre- 
senting the need for an interpolation in order to extract the result of the form factor, /^^, at q^ = 0. 

In order to reach small momentum transfers for the pion form factor and q^ = for the K ^ n 
form factors, we use partially twisted boundary conditions [11, 12], combining gauge field config- 
urations generated with sea quarks obeying periodic boundary conditions with valence quarks with 
twisted boundary conditions [11-17]. The valence quarks, q, satisfy 

qixk+L)=e"^'qixk), (^=1,2,3), (3.3) 

where 6 is the twisting angle. 




^3 
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Our method is decribed in detail in [6, 7] and proceeds by setting = for the spectator quark, 
denoted by ^3 in the above diagram. We are then able to vary the twisting angles, 0, and df, of the 
quarks before (^2) and after (^1) the insertion of the current, respectively. The momentum transfer 
between the initial and final state mesons is now 

<?2 = {Ei{piA)-Ef{pf, Bf)f - {{pi + Bi/L) - {pf + df/L)f , (3.4) 



where E{p, 0) = y + {p + d/L)^. Hence it is possible to choose 0,- and df such that = 0, 
which from now on we refer to as 6^ and 6,^ for when we twist a quark in the Kaon and Pion, 
respectively. 

In order to extract the matrix elements (3J_) from a lattice simulation, we consider ratios of 
three- and two-point correlation functions. For the pion form factor, we consider the ratios given 
in Eqs. (3.4) and (3.5) in [7], while for the K ^ n form factors, we consider the following ratios 



Ri,P,P,{pi,Pf) = 4,/EiEf^l hi,(t^,p^c,M^,pt) 



Cp,Pf {t,Pi,Pf) CpfP- {t,Pf,Pi) 



n ^ \ _ A / 77 17 Cp.p^{t,pi.pf) I Cp. {t,,nk-t:Pi)Cp^ {t.Pf)Cpf (4ink,P/) 

R3,P,Pf[PuPf) - ^^/Eitf c,^(,^_^^^^^,) y C.^,(,,,„,-r,p,.)Cp,(f,P,)Cp,(4.„,,p,) 



(3.5) 



We deviate slightly from the method outlined in [6] for extracting /^;j(0) from the ratios. 
Previously we considered only the time-component of the vector current and solved for f^j^{0) = 
fxK (0) linear combination 

/L(0) = ^"■f -f-) -R<^^Kn{Ap^{EK-m ) (^^ 2,3). (3.6) 

) {EK + m^){mK-E^)-{mK + E^){EK-m^) ^ ' ' ^ 

This, however, is just one of many expressions that can be obtained when we solve the system of 
simultaneous equations that are obtained when we consider all components of the vector current, 
Vu , rather than just V4 that was considered in [6] 

Ra,Kn0K,'^,VA) = i^K + m^) + f^^{Q) {Ek - m^) 

RaM^Ay^) = fiM{mK + En)+fKM{mK-E„) 

RaM^^ ^n.Vi) = fLiO) dni - IkM en.i • (3.7) 
We can now proceed to solve this overdetermined system of equations via minimisation. 



4. Pion form factor results 

In Fig. ^ we show our results for the form factor f^^{q^) for a pion with niji = 330 MeV for 
a range of values of both using periodic he's and partially twisted he's (set A and sets B&C 
respectively in the left plot of figure). The vertical dashed line indicates the smallest momentum 
transfer available on this lattice with periodic he's. The (blue) dashed line is the result of a pole- 
dominance fit to our data points, while the (red) dot-dashed curve is obtained from the result of 
QCDSF [18] evaluated at = 330 MeV. 
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Figure 2: f^^{q^) from a 24^ x 64 lattice with nin = 330 MeV using partially twisted bc's. 



On the right of Fig. ^ we have a zoom into the low = —q^ region. The triangles are our 
lattice data points for a pion with niji = 3 30 MeV, and the magenta diamonds are experimental data 
points for the physical pion. 

Because our values of are very small, we apply NLO chiral perturbation theory (ChPT). 
In NLO ChPT, the pion form factor depends only on a single low energy constant (LEC) (Lg for 
SU(3), or II for SU(2)) 

/su(2),NLo('7') = l + ^[-2llq^+4^^{ml,q\ll')] (4.1) 

/su(3)NLo(^') = 1 + 4 [4L5^2^4jr(m2,g2^At2)+2jf (4-2) 
h 



where 



and 



^im" q' u') = _ log !!!! (4 3) 



3^18 6 V X ^y^(^x-4)/x-\J 

with H(x) = —x/6 + 0{x^l^) for small x. Provided our pion mass is light enough, we can use the 
q^ dependence of f^^{q^) to extract this LEC. The grey dashed curve on the right hand of Fig. ^ 
shows our SU(2) fit to the nij^ = 330 MeV pion form factor data. 



Once the LEC is determined from this fit, we insert the physical pion mass in (^.1) to obtain 
the solid blue curve. In addition we also represent the PDG world average [2] for the charge radius 
using the black dashed line. Our best estimate for the pion charge radius comes from the SU(2) 
NLO ChPT fit to the three lowest points and is 

(r^) =0.418(31) fm^ . (4.5) 

The fact that our result is in agreement with experiment, (r|) = 0.452(II)fm^ [2], gives us confi- 
dence that we are in a regime where chiral perturbation theory is applicable. 
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5. Kit, form factor results 



As explained in Sec. ^ we calculate the K ^ n form factor directly at = by setting the 
Kaon and Pion in turn to be at rest, while twisting the other one such that = 0. We refer to these 
twist angles as dji and Ok, respectively. We then get the following equations: 



{KipK)\V^\7l{0)) 
{KiO)\V^,\7c{p^)) 

By simply solving the simultaneous equations 
for each of the jJ. components separately we 
find that the errors in /^;j(0) and f^j^{0), are 
much larger than the errors in the matrix ele- 
ments. We have managed to circumvent this 
by looking at all the /i components simultane- 
ously, and then performing a minimisation 
on the overdetermined system of equations to 
find the values of /^^(O) and f^ni^) that best 
fit the equations. 



fKKi^)PK,ti-fK7li^)PK„ 
fiMPn,ii+fKMPn,, 



(5.1) 



LIO 



To obtain the matrix elements (5.1), we 



consider different combinations of Ri and R^, 
(p3). We find that all combinations lead to 




Figure 3: K^t, form factor, f^^{q^), evaluated at 



q —0 directly using twisted boundary conditions. 
Results are compared with data at ^ and fits 
from [3] 



consistent results, with the best combination 
being that we use Rt, for all matrix elements 
except for the case where the pion is twisted 
and we are considering the 4* component of 

the vector current. Using this set up, we obtain our preliminary results for f^jiiO) and f^j^{0) (for 
a pion mass of nij^ = 330 MeV) 



/+,(0) =0.9742(41), /^,(0) = -0.113(12) . 



(5.2) 



Our result for f^j^{0) = /^^^(O) is shown in Fig. |3| where we compare with the previous determi- 
nations in [3] which used pole /^,^(0) = 0.9774(35) and quadratic /+„,/(0) = 0.9749(59) func- 
tions to interpolate between q^„^,^ and negative values of q^. In our previous result, f^jii'^) — 
0.9644(33)(34)(14), these were combined, taking a systematic error of (34) for the model depen- 
dence. This contribution to the error has been eliminated in our new calculation. 

We conclude that using partially twisted bc's for the Kej, form factor, is an improvement on 
the conventional method as it removes a source of systematic error, while keeping comparable 
statistical errors. Another source of systematic error in our result in [3] is due to the slight difference 
between our simulated strange quark mass (am^ + arrives — 0.043) and the physical strange quark 
{anis + am^sfi — 0.037) [8], and we are in the process of determining the effect this has on our result 
through a simulation with a partially quenched strange quark mass of anis + amres — 0.033. We 
also plan to combine our results with the latest expressions from chiral perturbation theory [19]. 
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